Super-Diffusive Spin- Transport as a Mechanism of Ultrafast Demagnetization 
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We propose a semi-classical model for femtosecond-laser induced demagnetization due to spin- 
polarized excited electron diffusion in the super-diffusive regime. Our approach treats the finite 
elapsed time and transport in space between multiple electronic collisions exactly, as well as the 
presence of several metal films in the sample. Solving the derived transport equation numerically 
we show that this mechanism accounts for the experimentally observed demagnetization within 
200 fs in Ni, without the need to invoke any angular momentum dissipation channel. 



Excitation with femtosecond laser pulses is known for 
more than a decade to cause an ultrafast quenching of the 
magnetization in metallic ferromagnets PQ . The achieved 
demagnetization times are typically 100-300 fs for ferro- 
magnets such as Ni [TJI2]- Hence, laser-induced demag- 
netization opens up new, interesting routes for magnetic 
recording with hitherto unprecedented speeds [3j. How- 
ever, in spite of the technological importance the mech- 
anism underlying the femtosecond demagnetization re- 
mains highly controversial. A common belief is that there 
should exist an ultrafast channel for the dissipation of 
spin angular momentum [JHH]- Several such mechanisms 
through which an excited electron can undergo a spin- 
flip in a ferromagnetic metal are currently being debated. 
The main proposed mechanisms for a fast spin-flip pro- 
cess are a Stoner excitation, an inelastic magnon scat- 
tering, an Elliott- Yafet-type of phonon scattering |H[S], 
spin- flip Coulomb scattering [6], laser- induced spin-flips 
[7], or relativistic quantum electrodynamic processes [5]. 
An effect that, until recently [9], has been regarded to 
play only a marginal role is the spin-polarized transport 
of laser-excited hot electrons. 

In this Letter we show that spin-dependent transport 
of laser-excited electrons provides a considerable con- 
tribution to the ultrafast demagnetization process and 
can even completely explain it. We demonstrate this by 
developing a transport equation for the super-diffusive 
flow of spin-polarized electrons. A few approaches to 
describe the electron motion have been attempted pre- 
viously [101 111] , In our theory, however, we take into 
account the whole process of multiple, spin-conserving 
electron scattering events and electron cascades created 
by inelastic electron scattering. Also the presence of dif- 
ferent metallic films in the probed material is treated. 
We solve the developed theory numerically for ferromag- 
netic Ni, for which the femtosecond demagnetization is 
well documented [TJ ® [T2] , and show that a large demag- 
netization in a few hundred femtoseconds is generated. 

The typical geometry for a femtosecond laser experi- 
ment is depicted in Fig. [I] The intense laser beam creates 
excited hot electrons in the ferromagnetic film, which 
will start to move in a random direction. Our goal is 
to compute the time-dependent magnetization resulting 
from the super-diffusive motion in the laser spot. Due 




FIG. 1: (Color online) Sketch of the super-diffusive processes 
caused by laser excitation. Different mean-free-paths for ma- 
jority and minority spin carriers are shown and also the gen- 
eration of a cascade of electrons after an inelastic scattering 
(elastic scatterings are not shown for simplicity). The inset 
shows the geometry for the calculation of the electron flux 
term in the continuity equation. 



to the fact that the electronic mean-free-path (up to a 
few tens of nm) is much smaller than the diameter of 
the laser spot normally used, see e.g. [2j, we developed a 
uniaxial model where only the z dependence is kept, but 
the formalism can easily be applied to describe systems 
with less symmetry. 

After the absorption of a photon, an electron will be 
excited typically from a d-band to the sp-likc bands above 
the Fermi level. The mobility of sp-like electrons is much 
larger than that of d electrons. Therefore in the follow- 
ing we treat the d electrons as quasi-localized and com- 
pute the transport only for the mobile sp-like electrons. 
We furthermore use that the optical excitation is spin 
conserving. Because of the very small linear momentum 
carried by a photon the angular probability density of 
the emission direction can be safely considered isotropic 
over all solid angles. Hence, the excited electron will 
start moving in a random direction. The outgoing trajec- 
tory is treated as a straight line up to the first scattering 
event. As the derivation involves many variables, in the 
following we show explicitly only those that are relevant 
at that moment. We start with describing the motion 
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of electrons before the first scattering event, which (for 
simplicity) we call first-generation electrons. The trajec- 
tory s(t) for a first-generation electron is given implicitly 
by f£ ds'/v (z (s')) = t, where v (a, E, z) is the velocity, 
z (s) the z-coordinate of the particle when the coordinate 
on the trajectory is s, and t the time needed to reach s. 
The probability to reach a point s without being scat- 
tered is P (s) = exp [— Jq ds'/r (z (s' j) v (z (a'))] , where 
r (a, E, z) is the lifetime. The velocity and lifetime are 
both considered dependent on the position z (therefore 
on the material) and on the spin a and energy E but not 
on the direction of motion. 

Using that the probability of the emission direction is 
isotropic and integrating over all possible emission angles, 
we compute the statistically averaged flux <f> at a time t 
through a surface located at z due to an electron that 
starts its movement at zq at time to (see inset in Fig. [I]), 
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and used that an electron which suffers a scattering does 
not give any contribution to the first-generation flux. 

If, instead of exciting a single electron, a distributed 
source of excited electrons is present, the total first- 
generation flux due to all electrons with spin a and energy 
E is 

/-foo rt 
dz / dt S ext (z ,t a )(f>{z,t;z ,t ), (3) 
-co J — oo 

where S ext — S ext (<r,E,z,t) is the electron source term 
that has to be computed from the spatial and temporal 
profile of the laser and the absorption probability. We 
define the operator 0, by 4>S ext = <E>. 

Once we have the expression for the flux it is straight- 
forward to write the continuity equation for the den- 
sity of first-generation electrons of given spin and energy, 
nW (a,E,z,t), 
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The second term is the scattering term that acts as a 
reaction term removing electrons from the density when 
they are scattered. 

To describe the motion of electrons after the first scat- 
tering event, we make the assumption that the angular 
probability density of the emission is again isotropic and 
uncorrelated to the incoming direction. In the case of 
scattering with phonons, impurities, and other large mass 



particles this is almost exactly true. Also for scattering 
of sp-like electrons with the much heavier d-electrons it 
is a suitable approximation since it leads to an underes- 
timation of the diffusive process. The density of second- 

§ieration electrons, ni 2 J (a,E,z,t), is described by Eq. 
, too, when instead of using S ext we use (a, E, z, t), 
defined as: 



,J T(a',E',z) 
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which is the scattering term coming from the first- 
generation weighted by the transition probability after 
a scattering, p (a, a' , E, E' ,z). The latter includes the ef- 
fect of both elastic and inelastic scatterings as well as 
the generation of cascade electrons. We note that in 
principle the transition probability could handle spin-flip 
events but we assume them to be negligible. We de- 
fine the operator S as Sn^ = S^. Applying the same 
procedure we obtain the equation for the density of the 
third-generation electrons rJ 3 l, and so on. Summing up 
everything we derive a set of coupled transport equations, 
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/ is the identity operator. 

Eq. describes the fast transport of laser-excited 
electrons. A major question is whether it can explain the 
femtosecond demagnetization observed in pump-probe 
magneto-optical experiments. Before addressing this 
question we stress that the developed transport process 
is both different from ballistic and diffusive transport. 
Standard diffusive processes that are governed by Brow- 
nian motion are characterized by the variance of the dis- 
placement of a particle distribution a 2 which grows lin- 
early with time: a 1 (t) = t 1 , with 7=1 pj . Ballistic 
diffusion is characterized by 7=2. The here- developed 
electron motion description is in the category of super- 
diffusive processes (7>1) with the further distinction that 
7 is time-dependent and goes from a ballistic regime 
7=2 for small times to normal diffusion 7=1 for long 
times. This furthermore emphasizes that a standard dif- 
fusion model is inapplicable to electron motion on the fs 
timescale. 

Super-diffusive transport may give rise to demagneti- 
zation, because, first, laser-excited electrons in sp-bands 
have high velocities (about 1 nm/fs), and second, ex- 
cited spin majority and minority electrons have different 
lifetimes. The latter lends excited majority carriers in 
typical 3c? ferromagnets a high mean-free-path whereas 
minority carriers are much less mobile. This may lead 
to a depletion of majority carriers in the magnetic film 
and a transfer of magnetization away from the surface. 
Moreover, an excited electron experiencing an inelastic 
scattering with another electron will transfer part of its 
energy to the other one, generating an electron cascade. 
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FIG. 2: (Color online) Calculated spatial magnetization pro- 
file of Ni at three times caused by laser excitation (at t=0 fs). 
The resulting magnetization profile is given by the full curve, 
the initial one by the dotted curve. The electric field correc- 
tion is illustrated by the dashed curve. The surface of the film 
with the vacuum is at nm depth, the Ni film extends up to 
15 nm depth, the remaining is the Al film. 



The newly excited electron will, now, have enough en- 
ergy to contribute with its motion to the demagnetiza- 
tion. The occurring energy transfer is computed here 
from the classical treatment of a two-particle collision. 

To assess how much demagnetization such process can 
generate we solve Eq. ^ numerically for the case of a 15- 
nm Ni film on an Al substrate. This system was recently 
investigated by Stamm et al. [5] . The energy dependence 
of the spin-polarized densities in Eq. Q is discretized in 
12 channels from up to 1.5 eV, the pump-laser energy. 
The interface with the vacuum is treated as a reflecting 
surface, but the excited carriers in the Ni film can pen- 
etrate the Al. A pump laser with a temporal profile of 
60 fs (full-width at half-maximum) is adopted. Further- 
more, excitation rates are taken from Ref. [2], the ratio 
of excited spin-up to spin-down electrons from |14j . and 
the spin-dependent inelastic lifetimes and velocities from 
Ref. [IS]. Unfortunately, Ref. jTHj does not provide in- 
elastic lifetimes for excited electrons with low energies 
(i.e., from the Fermi level to 0.5 eV above). We made a 
linear extrapolation to estimate lifetimes for the lowest 
energies; the influence of this approximation is discussed 
further below, together with those of other approxima- 
tions. The clastic lifetime contribution can be sample 
dependent. We used as elastic inverse lifetimes 100% of 
the inelastic ones. Values from 20% to 200% have given 
no relevant changes in the demagnetization. 

Fig.[2]presents computed spatial (z coordinate) profiles 
at three times. The maximum of the excitation laser is 
at fs, but due to its temporal profile (cf. Fig. [3]) there 
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FIG. 3: (Color online) Computed laser-induced demagnetiza- 
tion in Ni. The shaded area shows where the theoretical result 
is expected to be (depending on the inelastic lifetime). For 
comparison we also show recent experimental XMCD data [2]. 
The used time-structure of the laser pulse (in a.u.) is depicted 
by the red (online) line. 



is already an effect visible at t=0 fs. The laser excita- 
tion is spin-conserving, hence in itself it does not cause 
any magnetization change. But due to spin-dependent 
fast transport of excited electrons the magnetization (in 
fiB I atom) begins to deviate from the initial one in the 
surface region. This process continues to £=90 fs, where 
the magnetization in the near-surface area has become 
considerably reduced; the super-diffusive flow of spin-up 
electrons into the Al film causes a magnetization of the Al 
film. At £=300 fs, well after the pump-laser has vanished, 
the continuing motion of excited electrons has created a 
further demagnetization, but the densities of excited car- 
riers are already quite reduced. The latter stems from 
inelastic collisions which an excited electron suffers until 
it has lost its energy, i.e. the process of thermalization. 
Due to scattering in the Al layer, there is a back-flow 
from the Al that creates a minority spin accumulation 
on the Ni side, leading to a stronger demagnetization at 
the interface, visible for t=300 fs. 

The carrier flow creates charged regions that gener- 
ate an electric field. This field is not negligible for laser 
fluences that give demagnetizations of the order of tens 
of percent. It will, however, act equally on both spin 
channels and on all conduction electrons, not only the 
excited ones that are a small percentage. As an extreme 
scenario we computed the back-flow necessary to com- 
pensate the free charge generated by the super-diffusive 
transport; the results are shown in Fig. [2] Its influence 
on the magnetization profile is indeed small. 

Next, we study the calculated time-evolution of the 
magnetization in Ni and compare to results obtained in 
a recent pump-probe experiment [2], in which the x-ray 
magnetic circular dichroism (XMCD) was used to probe 
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the magnetization. The time-evolution of the magneti- 
zation is shown in Fig. [3] The computed magnetization 
is convoluted with the probing beam's temporal profile 
(FWHM=100 fs) and averaged over the Ni-film thick- 
ness, because the XMCD signal is sensitive to the Ni only 
and almost uniformly probes the whole film thickness. 
The predicted demagnetization is in reasonable agree- 
ment with the experimental XMCD data. After 200 fs 
the magnetization is reduced to about 50%. The XMCD 
experiment seems to indicate a somewhat faster demag- 
netization, but we must be careful when directly com- 
paring the absolute values of measured and calculated 
magnetization. It has been shown that, while demag- 
netization is occurring, the XMCD signal also contains 
artifacts due to state-blocking effects |16j . The shaded 
area gives the temporal behavior computed for two dif- 
ferent zero-energy extrapolations of the lifetime data |15j . 
A linear extrapolation of lifetimes for energies near the 
Fermi level gives an inelastic majority spin lifetime of 140 
fs; a rigid cut-off of the low-energy lifetime gives a value 
of 50 fs. The former corresponds to a stronger, the lat- 
ter to a weaker demagnetization (~40%). Using longer 
inelastic lifetimes does not lead to a notable increase of 
the demagnetization. 

We note that our model does not include two effects 
that might give additional contributions: the generation 
of Auger electrons and the response of the material to 
the new electronic state. The first effect provides a de- 
layed source of fast electrons, improving the effectiveness 
of the transport-driven demagnetization. The second ef- 
fect is more complicated. The inelastic lifetimes might 
be altered by the new electronic distribution which may 
reduce the effectiveness of the process. 

A confirmation of our theory is that it explains the 
striking difference between magnetic metals and mag- 
netic dielectrics. The latter exhibit much slower demag- 
netization times jTTl[T8] . As the super-diffusive transport 
processes are inhibited in insulators, the demagnetization 
has to evolve through the much slower dynamics of spin- 
flip processes. 

In contrast to Ni, a very slow laser-induced demagne- 
tization has been observed for metallic 4/ ferromagnets 
as Gd [TO] . Gd is different from metallic 3d ferromagnets 
because the spin-polarized 4/ states-which contribute 7 
/is to the Gd moment-are well localized and not reached 
by the pump laser. Only the Gd 5d band electrons near 
the Fermi energy are excited by the pump laser, but these 
electrons contribute only 0.55 fiB to the atomic moment. 
Our theory is consistent with the two-time-scales demag- 
netization in Gd, because the super-diffusive transport 
of laser-excited 5d electrons would create a small, fast 
5d demagnetization, but then the laser-imparted energy 
needs to be transferred from the d electrons to the lattice 
and the / electrons to achieve sizable demagnetization. 
Especially the latter process is slow (~80 ps, [EI]), as it 
involves spin-flip phonon scattering. 



Malinowski et al. [5] observed a demagnetization con- 
tribution due to spin transport, however, they did not 
have a detailed description of the process to quantify its 
amount. We have developed here a model that can ex- 
plain laser-driven fs demagnetization on the basis of spin 
angular-momentum conserving super-diffusive transport. 
Numerical solutions of the transport equations show that 
a substantial demagnetization of ~50% is created within 
200 fs, without the need to invoke any spin-flip chan- 
nels. At this stage we cannot exclude other demagneti- 
zation mechanisms, but our calculations show that super- 
diffusive processes play a main role, and can explain 
wholly the ultrafast demagnetization process during the 
first few hundred femtoseconds. On longer timescales 
other well-known effects that are not included in this 
treatment but that can be modeled by Landau-Lifshitz- 
Gilbert dynamics |20j come into play. A combination 
of these approaches is required to treat the full fs to ps 
time-domain. A test of our theory could be to observe 
the flow of spin-polarized electrons in the non-magnetic 
substrate, which, for sufficiently thin substrates could be 
detectable by magneto-optical techniques. 
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